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Abstract 

We present a new algorithms to discretize a decoupled forward backward stochastic 
differential equations driven by pure jump Levy process (FBSDEL in short). The 
method is built in two steps. Firstly, we approximate the FBSDEL by a forward 
backward stochastic differential equations driven by a Brownian motion and Poisson 
process (FBSDEBP in short), in which we replace the small jumps by a Brownian 
motion. Then, we prove the convergence of the approximation when the size of small 
jumps £ goes to 0. In the second step, we obtain the L'p Holder continuity of the solution 
of FBSDEBP and we construct two numerical schemes for this FBSDEBP. Based on 
the L'P Holder estimate, we prove the convergence of the scheme when the number of 
time steps n goes to infinity. Combining these two steps leads to prove the convergence 
of numerical schemes to the solution of FBSDEL. 
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1 Introduction and summary 



In this paper, we are concerned by discretization of a system of decoupled forward-backward 
stochastic differential equation (FBSDEs in short) driven by a pure jump Levy process 

= Xo + /o*6(X,)dr + /o7j,/3(X,-)M(de,dr), 
= g{XT) + f f{Qr)dr - VrM{de, dr) ^ ' ^ 

Here 6 := (^X^Y, p{e)Vei'{de)) and M{E,t) = J^^j^ e/x(de, lir) where fl{de,dr) := 
fi{de, dr) — ^[dejdr an independent compensated Poisson measure and a Poisson random 
measure on M x [0, T] with intensity v satisfying / I A |epz^((ie) < oo. 

Numerical discretization schemes for FBSDE have been studied by many authors. In the 
no-jump case, Ma et al. [21] developed the first step algorithm to solve a class of general 
forward-backward SDE. Douglas et al. [14j suggest a finite difference approximation of the 
associated PDE. Other discrete scheme have been considered in [7], [S] and [llj mainly 
based on approximation of the Brownian motion by some discrete process, Gobet et al. 
|18j proposed an adapted Longstaff and Schwartz algorithm based on non-parametric re- 
gressions. In the jump case, to our knowledge, there is only the work of Bouchard and Elie 
[5] in which the authors propose a Monte-Carlo methods in the case when v{U.) < oo. 

The main motivation to study the numerical scheme of a systems of above form, is to treat 
the case when z^(M) = oo, which means the existence of an infinite number of jumps in every 
interval of non-zero length a.s.. In this sense, we should mention the important work on 
the approximation of stochastic differential equation studied by Kohatsu-Higa and Tankov 

m. 



Since we are interested in the case of z^(M) = oo, we will follow the idea of |20j to ap- 
proximate (jl.2p without cutoff the small jumps smaller than e, which should improve the 
approximation scheme. Then by using the approximation result of Asmussen and Rosinski 
[2] we replace the small jumps of the driven-Levy process with a{e)W where is a stan- 
dard Brownian motion and (T^(e) := e'^v{de). 

In the aim to approximate (|l.ip . we cut the jumps at e as the following 



Xt = Xq + f^b{Xr)dr + jlp{X,-)dRr + fj'j, P{X,-)M{de,dr) 



Yt = g{XT) + !^ f{er)dr - // VrdRr - // Je VrM{de, dr) 



(1.2) 



where Rt = f^^^^^eM{de,dr), E" := {e G M, s.t / |e| < e}, Ee := {e G M, s.t / |e| > e} 
and E :=R = E^U Ee- 

The idea we propose is to discretize the solution of (jl.ip in two steps. In the first step, we 
approximate (jl.2p by the following FBSDE: 



= X^,+JibiX^,)dr + j;^/3{X^Me)dWr + j;,fj,JiX^.^)M{de,dr) 
Yf = g{X^) + f{@^,)dr - Z^JWs - L U^{e)M{de, dr) 
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Here 0^ := (^X^,y^,r^j and := p{e)U'^{e)ev{de). Further, we show that for a finite 
measure m defined by m{E) := f^e'^i'{de), our error 



Err^(Y,V) 













2' 


E 


supiy* -y/p 


+ E 


sup 


[ VrdRr- [ Z^dWr 






t<T 




t<T 


Jo Jo 





+E 



T 

JE^ 



\Vr - U^.ie)\'^m{de)dr 



is controlled by o"(e)^, which means that the solution of (I1.3P converges to the solution of 
(jl.ip . as the size of small jumps e goes to (See Remark 12. ip . We also derive the upper 
bound 



E 



sup \Xt - XI 

t<T 



ei2 



(1.4) 



The second step consists of discretizating the approximated FBSDE (jl.Sp and studying its 
convergence to (II. 2p . For this purpose we consider two numerical schemes, the first one is 
based on discrete-time approximation of decoupled FBSDE derived by Bouchard and Elie 
[5]. More precisely, for a fixed e, given a regular grid vr = {tj = iT/n,i = 0, l,...,n.}, 
the authors approximate by its Euler scheme X'^ and {Y^ , ,T^) by the discrete-time 
process {Yf^,Z^,T^) 



= Xl + ^b{Xrj + P{X[Jaie)AW,+, + J^JiXrjM{de,{t^ 

= nEfe^At^,+i/7-i,] 

= nE [y,-^^ Jj,^ p{e)M{de, {U,U+i])/Tt, 

= E\Y^JT^+y{xl,^7:,W 



VTT 



(1.5) 



on each interval [tj,tj+i), where the terminal value := g{X[^). Under Lipschitz conti- 
nuity of the solution, the authors proved that the discretization error 



Err,(Y',Z',T' 



supE 

t<T 



\Yf-Y,^\']+ rE[|Zf-ZrP + |Ff-f-p]dt (1.6) 
Jo 



achieves the optimal convergence rate n Finally, we derive the first main result of this 
paper in Proposition 13.11 showing that the approximation-discretization error 



Err^^^,){Y,V) := supE 



t<T 



\Yt - y/P + 



VfdRr 



* _ 2 
Z^dWr 



+ lir 



=177112 
11^2 



(1-7) 



is bounded by C{n ^ + cr(e)^) and converges to as (n,e) tends to (oo,0), where F := 
p[e)Veu{de) . Taking e = n"^/^, our approximation-discretization achieves the optimal 



convergence rate n 



-1/2 



The second numerical scheme has been inspired from the paper of Hu, Nualart and Song 
[19j . Where the authors study a backward stochastic differential equation driven by a 
Brownian motion with general terminal variable ^. They propose a new scheme using the 
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representation of as the trace of the Malhavin derivatives of Y^. Their discretization 
scheme is based on the L^-Holder continuity of the solution Z^, to obtain an estimate of 
the form 

E|Zf <K\t-s\^, 

which imphes the existence of a 7-H61der continuous version of the process for any 
7 < ^ — i. In this sense, our article extend the work done in [19] to a forward-backward 
stochastic differential equation with jumps and terminal value g{Xj,). Similarly to [19], we 
obtain the following regularity of F*" 

E|rf-n|p < c|t-s|^, 

which allows us to deduce the existence of a 7-Holder continuous version of the process 
r*" for any 7 < ^ — ^- Finally, on one hand, we use the representation of Z^ and F^ as 
the trace of Malliavin derivative of to derive our a new extended discretization scheme 
for the solution (Y^,Z'^,T^) of (jl.3p . From other hand we approximate by X'^ the 
continuous-time version of the Euler scheme, that is for a fixed e > 

' X[ = X;^+aie)b{X;^){t-cPf) + a{e)P{X;„){Wt-W^^) + J^J{X;^)M{de,{^^ 

= E£[^^^,jMx^)Dt,x^ + ElZl£r.,.,t,^Af(®^^^^^^ 

(1.8) 

with terminal values Y,l = g{X^), Z^ = a{e)d,g{X^)^{X^) and Ul^, = g{X^ + ^{X^)) - 
g{X^), where £ti,tj ^'^d Stl^j detailed in section 4. 

The key-ingredient for computation of discretization error, is based on the L^-Holder con- 
tinuity of the solution (y^, Z^,T^). This allows us to prove that 



Errl{Y^,Z'',T'''^ 



E max 

0<i<n 



Yf - y.'^P + \Zl - ZT\^ + |F^ - Ff 



TT |2 I \T^S l~i7r |2 



is controlled by |7r| ^''^M . Then we obtain the second main result of this article in Theorem 



131 which stating that 



Errl,{Y,V) := max sup [\Yt - Y^?] + 



/ VrdRr-Y^Zl^Wt 



is of the order o"(e)^ + |7r| '"""R and converges to as the discretization step (e,n) tends 
to (0, 00). 



4 



The importance of the above scheme, is it can be adapted to the case when the a terminal 
value is not given by the forward diffusion equation , as it 's the case in |19] . However, 
this scheme remains to be further investigated. 

The two numerical schemes above are not directly implemented in practice and require an 
important procedure to simulate the conditional expectation. However, there exist different 
technics which can be adapted to our setting to compute this conditional expectation and 
we shall only mention the papers: [3], [6] , [9] and [18]. 

The paper is organized as follows. In Section 2, we prove the convergence of the approxi- 
mated scheme. In Section 3, we describe discrete-time scheme introduced in [5] and state 
our first main convergence result. In section 4, we extend the new discrete scheme of [H] 
and state our second main result. We also discuss a general case of BSDE. Section 5, is 
devoted to Malliavin calculus for a class of FBSDE with jumps, we then get the L^-Holder 
continuity of and via the trace of the Malliavin derivatives of . 



2 Approximation of decoupled FBSDE driven by pure jump 
Levy processes 

Let (ri, J-", F = (J-')t<7^, P) be a stochastic basis such that J-q contains the P-null sets, J-t = J' 
and F satisfies the usual assumptions. We assume that F is generated by a one-dimensional 
Brownian motion W and an independent Poisson measure ^ on [0, T] x E. We denote 
by F^ = {J-t^)t<T (resp. F^ = {J't)t<T) the P-augmentation of the natural filtration 
of W (resp. n). As usual, we denote by B{X) the Borel set of topological set X. We 
introduce the following subset: E^ := {e G M, s.t / |e| < e}, E^ := {e G M, s.t / |e| > e}, 
E:=R = E^\J Ee. 

The martingale measure /x is the compensated measure corresponding to Poisson random 
measure /x, such that Ji{de, dr) = fx{de, dr) — v{de)dr, where is a Levy measure on E 
endowed with its Borel tribe £. The Levy measure v will be assumed to satisfy i/(R) = oo 
and J]g |epz^(de) < oo. Throughout this paper we deal with the measure M defined by 

M{t,B)= [ efL{dr,de), B e B{E) 

J[0,t]xB 

which can be considered as a compensated Poisson random measure on [0, T] x E and 
JjQ jj^^ e//((ir, de) is a compound Poisson random variable. We associate to M the cr-finite 
measure 

m{B) := / e^iy{de) B G 13(E). (2.1) 
Jb 

In particular, we have cr{e)'^ = m{E^). 

The measure M is taken to drive the jump noise instead of /i, in the aim to adopt the 
concept of Malliavin calculus on the canonical Levy space from [12] . 

For some constant K > 0, we consider four i^-Lipschitz functions with bounded derivatives 
/3 : R ^ M, 6 : M ^ M , 5 : M ^ M and / : Q X M X M X L2(£', I/, R) ^ R, where the first 
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derivative of b, /3 and g are bounded. 

Define p to be a measurable function p :£'—)■ M such that: 

sup \p{e)\ < K. 

For any p >2 we consider the following class of processes: 

• is the set of real valued adapted rcll process Y such that: 

||F||5P :=E ( sup \Yt\A < oo. 

\o<t<T J 

• is the set of progressively measurable M- valued processes Z such that: 



(2.2) 



\Z\ 



HP 



E 



iZJ^dr 



< oo. 



is the set of V £ measurable map U : Q x [0,T] x E M. such that: 

1 

rT 



\U\ 



LP 



E I / \Ur{e)\Pu{de)dr ] < oo. 
Je 



\Y\]P + -4- \]Tr\\p ^ " 

I WSP ' W'^W HP ' 1 1 1 1 LP 



The space := x x is endowed with the norm 
\\{Y,Z,U)\\t3P :-- 

M^'P the class of square integrable random variable F of the form: 

F = EF+ f UrdWr+ [ [ ij{r,e)fl{de,dr), 
Jo Jo Je 

where u (resp. ip) is a progressively measurable (resp. measurable) process satisfying 
sup«r E|iif I*' < oo (resp. sup^<2-E \ip{t,e)\Pi'{de) < oo). 



2.1 Approximation scheme 

In this subsection, we show that the approximation error 
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sup|yt-F/p 


+ E 


sup 


[ VrdRr- [ Z^.dWr 






t<T 




t<T 


Jo Jo 





+E 



T 

JEe 



Vr - U^{e)\'^m{de)dr 



converges to as e goes to 0. 



Theorem 2.1 Under the space {Q,,T,F), 

1. There exist a solution X on [0,T] of 

Xt = Xo+ [ b{Xr)dr+ [ P{X,-)dLr, 
Jo Jo 

where Xq G R. 

2. There exist a solution X^ on [0, T] of 



(2.3) 



= X'o+ [\{X^.)dr+ f p{X',)a{e)dWr+ f [ 
Jo Jo Jo Jei 



P{Xl.)M{dr,de),{2A) 



where Xq G 



Moreover 



E sup \Xt-Xl\'^ < Ca{ef. 
\o<t<r J 

For the proof, we state firstly the following Lemma 

Lemma 2.1 On the space (fi, J-", P), fixing e > 0, we have for p > 2: 



(2.5) 



E I sup \Xt\P < oo. 

,0<i<T / 



E sup \Xf\P < oo. 

\0<t<T I 



(2.6) 



(2.7) 



Proof. We denote by C a constant whose value may change from line to line. Using 
Jensen's inequality, Burkholder-Davis-Gundy inequality and Lipschitz property of b and (3 
we have: 



E sup \Xs\P < C sup ( E 

0<s<t 0<s<t 



< C \Xq\p+ / E 



Xo|P+(^y b{Xr)drj +J J (3{Xr)M{dr,de) 
t 







\b{Xo)\ + \Xr\ + \X^ 



dr 



+ f h 

10 Je 



|/3(Xo)| + |Xo| + |X^| ePu{de)dr 



< c ( |Xor + |6(Xo)r + \f3{Xo)T + E 



sup \Xu\Pdr 

0<u<r 



where C depends on t, b{XQ) and /3(Xo). We conclude the first assertion by Gronwall's 
Lemma. 



Following the same arguments, we obtain the second assertion. 



□ 



Proof of Theorem 12.11 The existence and uniqueness of such SDKs was already studied 
in the literature see e.g. [T7] and [1]. Then it remains to prove the estimate (j2.5p . 
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Using Jensen's inequality leads to : 



E sup |X„ -X^p < C E 

0<u<t \ 



+E 



\b{Xr) - b{X'^)\dr 







-1 2 



-E 



|/3(X,)-/3(X,^)|M(dr,de) 

Je. 

\l3{Xr)\M{de,dr)) + ( f \(3{X',)a{e)\dWr 



JE' 







By Burkholder-Davis-Gundy inequality, we get 



E sup \Xu-Xl\' < C(tE 

0<u<t 



{b{Xr)-b{X',)fdr 



+E 
+ E 

By Lipscliitz property of b and (3 



{p{Xr) - p{X^,.)f m{de)dr 



J Ee 

t 



l3{Xrfm{de)dr 



JE' 



+ CK 



PiXlfaiefdr 



E sup \Xu-Xl\^ < 

0<u<t 



< c 



J^E{Xr- X^f dr 

+<^f^f^ [/?'(^o) + /3'(^o') + + l^rf] dr 
dr + a{ef 



sup \6Xu\ 

u<r 



The result follows from Gronwall's Lemma. □ 
Finally, we can now state the main result of this section. 
Theorem 2.2 Under the space (J], J-", P), 

1. There exist a unique pair {Y,V) eS^ X H^, which solves the BSDE: 

Yt = giXT)+ r fi@r)dr- r VrdLr, (2.8) 
Jt Jt 

where 9 := (^X, Y, p{e)Veu{de)^ . 

2. For a fixed e > 0, There exist a unique solution (Y^ , Z'^ ,U^) € of the following 
BSDE: 

Yf = g{X^) + r f{Q'^)dr- r Z'^dWs - r [ U^{e)M{dr,de). (2.9) 

Jt Jt Jt J E^ 

with := (xsy^r^) andr' := p{e)U^{e)eu{de). 



Moreover, if supi<2^E|Vtp < oo, then there exist a constant C such that: 



Errl{Y,V) < Ca{ef. 



(2.10) 



Remark 2.1 Observe that 



E 



sup \Yt - Yi 

t<T 



e\2 



+ E 



sup 

t<T 



T 



VrM{de, dr) 



Z^dWr 



T 

JE^ 



W,{e)M{de,dr) 



< Errl{Y,V) 

Which shows clearly the convergence of the approximated scheme to lil.l]) . 

Proof, of Theorem 12.21 Existence and uniqueness of the solutions of BSDEs ()2.8p and 
(I2.9P was already proved, see e.g [4j. 

We are going to prove inequality (|2.10p . By Ito's formula applied to \5Y\'^ := \Y — Y^\'^ 
yields : 

\6Yt\'^ + / Zfdr+ I [ {U^.{e)-Vrfm{de)dr 

= |5(Xt) - 5(^t)|' + He? £ Vrdr - 5YrZ',dWr 

+2 r 5Yr if (9,) - / (9^)) dr + 2 r [ \{5Y,- + Vrf - Jy/J M(&, dr) 

Jt Jt J E^ L ^ 



T 



t J E, 



{5Y,- + U^^{e) - Vrf - 5Y^. M{de, dr). 



(2.11) 



Taking expectation in both hand-side of the above equality we get 



E 



T 



SY^+ I Zfdr + 



T 



t JE, 



{W,{e)-VrYm{de)dr 



E 



\g{XT) - g{X'T)? + a{e)^ r V^dr + 2 f 5Yr {f {Qr) - f m) dr 

Jt Jt 



From Lemma |2.H Lipschitz property of g and Jensen inequality we obtain: 



E 



Yt-YF\'+ Z'tdr + 



T 

t JEe 



{U^{e) -Vr)^m{de)dr 
< C (^a(e)2 + KeJ' (Yr - y/)(|X, - X'^\ + |y, - y/|)(ir 



+KK 
+ KE 



T ^ 



(Yr - y/) / p{e)e\U^.{e) - VMde) 



dr 



dr 



(Yr - y/) / p{e)Vreu{de) 



Using the fact that ab < aa^ + for some a > 0, yields to 



E 



\Yt - + /' Z'ldr + r [ iW,ie) - Vrf m{de)dr 
< C (^a{ef + K{l + a^ + + rf)¥. j\Yr-Y^^fdr + ^-K \Xr - X'^\^dr 
+ [ p{ef\U^(e)-Vr\'^mide)dr + ^E[ [ p{efV^^m{de)dr\, 

T Jt Je, T Jt Jo<\e\<e J 

where a and 7 are two constants taken such that ^ = = ^) we then get 

\Yt - y/|2 + r Z^ldr + r I {U',{e) - Vrf m{de)dr 
Jt Jt J E^ ^ 

< C (^{ef + {K + 2K^ + 2KK^)E {Yr - Y^fdr^ . (2.12) 



E 



Using Gronwall's Lemma, we deduce that 

E\Yt-Yl\^ < C(T{ef. 
Phigging this estimate in the previous upper bound, we get 

E / Zfdr + E [ [ {Uf.{e) - Vrf m{de)dr < Ca{ef 
Jo Jo Je^ 



(2.13) 



(2.14) 



Then 



E|5ytp + E /" Zfdr + E [ [ {U^{e) - Vrf m{de)dr < Ca{ef . (2.15) 
Jo Jo Je. 



Now using Burkholder-Davis-Gundy inequahty, we have 

rT f-T 



Esup|(5yj|2 + E / Zfdr + E [ [ {U^{e) -Vrf m{de)dr < Ca{ef . 

t<T Jo Jo Je, 

From other side, it follows by Burkholder-Davis-Gundy inequality and (|2.16p that: 



(2.16) 
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[ Z^.dWr 
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f VrdRr 






t<T 


Jo Jo 






t<T 


Jo 


t<T 


Jo 





< CE 



/ Zfdr + a{e^) [ Vfdr 
Jo Jo 

< E [ Zfdr + Ca{e^). 
Jo 

The result now follows by combining ()2.16p and (j2.17p . 



(2.17) 

□ 
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Remark 2.2 One can show that sup^<2- E| |^ < oo, following the same arguments as in 
Remark 2.8 in 119^ - 

Remark 2.3 In the BSDE i2.9\) . each time we change e, there exist a unique pair {Z^, U^) 
of predictable process, such that the BSDE i2. 9\) has a solution. 



3 Forward-backward Euler scheme 

In this section, we discretize the solution of ([L3]) by (X'^, ^ , Z'", T'^) de- 

fined by induction in (j3.3p and then we show the convergence of (X'^ ,Y'^ , Z'^ jV^) to the 
solution of (|1.2p . Thus let us recall some definition and notation. 
For each t G we define: 



and 



Zt = nE 



ZZ = nE 



ti+i 



Zsds/Tti 



nK 



ti+i 



Tsds/Tt^ 



rU+i 



Zlds/Tu 



IJt 



r rt 



rj =nE 



T:ds/Tu 



(3.1) 



(3.2) 



The process Zt^ and Ft; (resp. Z^. and TJJ can be interpreted as the best approximation 
of Zt^ and Fj, (resp. Z[^ and F^J. We know from Bouchard and Elie [5], that FBSDE (fZHjl 
has a backward Euler scheme taking the form: 



Z[ 

^ t 



XI + ^b{Xrj + a{e)AWi+, + /^^ P{Xl)M{de, (t^, t.+i]) 
nE fe^AW,+i/J-i,] 

nE [y-^^ /^^ p(e)M(de, (t^, ti+iD/J-*, 



(3.3) 



for which the discretization error: 



Errn{Y',Z',T') : = 



supEyy/-y,n'] + 11^' 

t<T 



" Z^||^2 + ||F^ -F^"2 



//2 



< Cn-1/2, 



(3.4) 



converges to as the discretization step ^ tends to 0. Means that the discretization scheme 
(13. 3p achieves the optimal convergence rate n"^/^. The regularity of Z^ and F*" has been 
studied in LP' sense in (5] when the terminal value is a functional of forward diffusion. 

It is well known also that 



maxE 



sup 

iS[ti,ii+i] 



7r|2 
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Our aim in this part, is to show that the approximation-discretization error between BSDE 
(fLljl and (fOl : 

2 



Err1^.{Y,V) := supE [jF^ - F,"|2] + supE 



(n,e) 



t<T 



t<T 







t 

VrdRr- I Z^dWr 




|r-f-||^2, 



(3.5) 



converges to as (e,n) — > (0, oo). 
The first main result of this paper is: 

Proposition 3.1 Assuming Lipschitz property of coefficients b and (3, the approximation- 
discretization error defined in /13. 5|) is bounded by: 



Err(^^,^,){Y,V)<c(n-'/^ + a{e)) . 



(3.6) 



Means that: 



Err^^,,){Y,V) — > 0. 

(n,£)--i>{oo,0) 



Proof. From ()3.5p . Jensen inequality and Burkholder-Davis-Gundy inequality, we have: 



Errl^^,){Y,V) < C (^snpE [\Yt - Yf \' + \Yf - Y,^ f] + \\T - T W^, + \\r' - 



■ t<T 



+ supE 

t<T 



Zf.dWr 



[ VrdRr - 

Jo 

From other side, it follows from Holder inequality that: 



p{e)e{Vr - U^.{e))u{de) <K{ j (K - U^{e)ym{de) ] v[E, 



(3.7) 



(3., 



Recalling that I'iE^) < cx) so that has a.s. only a finite number of big jumps on [0, T]. 
Combining the two last inequalities with (13. 4p leads to: 



Err^(^^)(Y,V) < C f + supE [ly* - F/p] + supE / 

\^ t<T t<T Jo 

+ [ [ E\Vr\'^m{de)dr + E 
Jo Je^ 



VrdRr - / ZrdW, 



T 



UO JE, 



\U^{e) - Vr\'^m{de)dr 



By Theorem 12.21 we get: 



ErrrnAY,V) < C (n^'^^ + aie)) , 



(3.9) 

□ 



where C depends on K. 

Remark 3.1 In the general case, as we neglect the small jump, the Brownian part in il.3\ 
disappears. In this case the assertion i3. 6\) can be replaced by: 

E^(^n,e){y,V)<Cn-^'\ 

Remark 3.2 Taking e = n~^/^, we obtain the optimal convergence rate n~^/^ in 

E^^r^,e){y,V)<Cn~^'\ 

which is exactly the approximation error in 
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4 A discrete scheme via Malliavin derivatives 



In this section, we generalize the new discrete scheme recently introduced by Hu, Nualard 
and Song [19] from a general BSDE, to our framework of decoupled forward-backward 
SDEs with jumps. For this aim, we use the Malliavin derivatives of Y to derive the discrete 
scheme. We first fix a regular grid vr := {ij := iT/n, i = 0, n} on [0, T] and approximate 
the forward SDE X'^ in (ll.3p by its Euler scheme X'^ already defined in ()3.3p . It'is hard 
to prove existence and convergence of Malliavin derivatives of X'^ . However, to avoid this 
problem, we can instead consider the continuous-time version of the Euler scheme, then we 
define the function (p for each t G [0,T]: 

:= max{tj, i = 0, n. / U < t}, (4.1) 

for which we associate: 

X[ := Xl^+a{e)b{Xl^){t-4><i) + a{e)P{Xl^){Wt-W^r.)+ f (3{X;„)M{de,{t,cP^]). 



It could be written as 

rt rt rt 



X^ := Xo+ [ b{X;r,)dr+ [ a{e)P{Xl^)dWr+ f [ /3(X;„)M(de, dr). 
Jo Jo Jo Je^ 



10 Jo Jo JE, 

It is well known that under Lipschitz property of the coefficients 

i/p 

sup \Xt-X^f 

iG[0,T] 



E 



The Malliavin derivatives of the continuous-time version of Euler scheme for 6 < s a.e. are: 
DeX[ = 1^ d,b{X;„)DeX;^dr + 1^ 1^ d^f3{X;^)DeX;„M{de,dr)+a{e)fi{X^ 

+aie) f d^P{Xl^)DeXl^dWr, 

De,eX: = [' De,eb{X:^)dr + [' [ De,eP{X;)M{de,dr) + a{e) f De^eli{X;)dWr 
Je Je Je^ Je 

We introduce some additional assumptions: 
(Al) f{t,y,^) doesn't depend on x. 

(A2) The first derivative of b and /? and g is a K-Lipschitz function 

\b'{x) - b\y)\ + |/3'(x) - /3'{y)\ + \g'{x) - g'{y)\ < K\x - y\. 

(A3) f{t,y,^) is linear with respect to t, y and 7. Moreover, there exist three bounded 
functions /i, /2 and /s such that : 

f{t,y,u) = hit) + /2(t)y + /3(t)7. (4.3) 
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Lemma 4.1 Under Lipschitz continuity of b and (3, we have for any q>l 

< oo, 



sup sup E 

0<e<Tn>l 



sup WDgX^fl 

e<t<T 



sup sup E 

0<e<Tn>l 



sup pe,e^rf' 
e<t<T 



< oo. 



(4.4) 
(4.5) 



For the proof see the Appendix. 

We then derive the fohowing theorem 

Theorem 4.1 Under assumption (A2), Lipschitz continuity of b and (3 and for any p > 2, 
we have, 



E 



sup \DgXf - DeX[\P 

tG[0,T] 



E 



sup \De,eX! - De,eX[\P 
te[o,T] 



i/p 



i/p 



(4.6) 
(4.7) 



Proof. Using Burkholder-Davis-Gundy inequahty, Jensen inequahty, inequahty (|4.2 
and Lemma |4. II 



E 



sup \DeXf-DeX[f 

sG[0,t] 







< CpE 


[/ 




Jo 



d.,b{XJr,)DeX;r. - d^b{Xr)DeX', 



dr 



+ 



JE. 



d,f3iX;„)DgX2n - d,f3iX',)DgX', v{de)dr 



2 \ 

dr I 



p/2 



which leads to 



< CpE 



\DeX, 



e\P 



\d^b{Xl) - d,b{Xl^)\P + W{X^,.) - d,p{xi. 



dr 



+ / \DeX',-DeXl 
Jo 

+ \x^-xi\p 



\dMx;.)\p + \d,Pix;.)\p 



dr 



< CpE (71^/2+ f sup \DeX^- DeXZl^dr 

\ Jo ue[0,r] ^ 



We conclude by using Gronwall's Lemma. Following the same arguments, we prove the 
second assertion. 

□ 



Now we derive the discrete scheme using the expression of and [/^ as the trace of 
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Malliavin derivatives of Y. From equation (jS.lSp and (j5.17p . the two Malliavin derivatives 
D^Yf , Dq^Yi could be expressed as: 

cT 
Jt 



(4.8) 
(4.9) 



where 

£t,r 



exp 



dyfi@l) 



d^f\Ql)p\e)m{de) du 



+ 




t JEs 



d^fiei)p{e)M{de,du) 



'-'t,r 



exp 



ae,u - ^ / p^{e)m{de) 



du + 




t je. 



Oie,uP{e)M{de, du) 



and 



Thus, we define our discrete scheme for i = n — 1, 1,0. and t £ [tj, tj+i) by induction 

Y,l^ + fiei^JAU^^/Ft^^ 

£r,^,,t,A9{xm.x^ + EkZl £i^,,,^,d.f{ei^,)Dt^xi^^/^tu/Fu 

(4.10) 

with terminal conditions 

Y,l = g{X^), Zl^ = a{e)dMX^)(3{X^), Ul^, = g{X^ + /3(Xf )) - ^(Xf ), 
where for any < i < j < n. 





= E 




= E 




= E 



1,1, 1, J 



exp 



£f 



\ k=i 

i-1 rt 

+ E 

k=i 

exp i 2^ 

I k=i 



9yf(®V-lf d,f{el)p^{e)m{de) 
^ Je^ 

d^f{Ql)p{e)M{de,dr)\, 



dr 



k+l 
tk J Eg 



\r,tk 



(4.11) 



^«''e,r,tfe / p\e)m{de) 



dr 



k=i ''^k 



tk + 1 



a0 j. -i.^p{e)M{de, dr) 



(4.12) 



and 



'^e,r,tk 



^tk 



f{ei+De,e@l)-f{&l) ^ 

X ir 



Ul^,p{e)u{de), 
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We are going to compute the discretization error of our discrete scheme and prove the 
convergence of the above scheme. We recall the expression of the error between the solution 
of ([OD and 



E max \\YF - y,^|P + iZf - ZJ\p + IH - 1^1 , 



where, = J^^ p{e)Ul^^u{de). 

We also recall the expression of discretization-approximation error between ()1.2p and (j4.10p 



Errl,{Y,V) := max sup E [jy^ - y^p] + E 



n-1 

+E 

i=0 



We conclude this section with the following Theorems whose proof are at the end of section 
5. 

Theorem 4.2 Under assumption \5. 11 we assume the existence of a constant L3 > such 
that : 



\f{t2,y,u) - f{ti,y,u)\ < L3\t2-ti\2. 

Then there exist a positive constant C independent of vr such that: 



ErrP{Y^, Z^ ^ ) < Cp|7r| '""^ . 

The second main result of this paper is summarized in the following theorem: 
Theorem 4.3 Under the same assumptions as Theorem \^.S\ we have 



(4.13) 



(4.14) 



ErrnAY.y)<C\ <y{e) + \n\ 



2log 



(4.15) 



Remark 4.1 The importance of the above scheme is it can be adapted to a backward SDE 
when the generator does not depends on the terminal value of a forward equation. 
Consider the following backward stochastic differential equation driven by pure jump Levy 
processes 

Yt = i+ f{r,Yr,j p{e)Vrei^{de)jdr- VrdLr (4.16) 
Which we approximate by 

y/ = i+ f f(r,Y^', [ p{e)Vreu{de))dr- [ Z'^dW^ - [ [ U^ie)M{dr,de) 

Jt ^ JEe ^ Jt Jt JEe 

(4.17) 
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We finally propose the below discrete time scheme, defined by terminal values Y^^ = ^, 
Zl = Dri and Ul = Dr^ei 



Zl 



(4.18) 



with ei = {r,Y^i,rj). 



Under the same assumptions of Theorem \4.^ we prove the convergence of the system 
to BSDE ^.16 ). Moreover, we obtain the upper bound 



max sup EHyt-y/Pl+E 



/ VrdRr-Y^ZlI^Wu +J2 I 
•'^ i=0 i=0 



E\Tt - ril'^dt 



<Cia\e) + \7T\ 



logj 



(4.19) 



5 Malliavin calculus for FBSDEs 

For ease of notations, we shall denote throughout this section the process {X^ ,Y^ , Z^ ,T^) 
by {X,Y,Z,T). 

In this section, we study some regularity properties of the solution {X,Y, Z,T). We recall 
the system ()1.3p using the new notations 

Xt = Xo + f;,biXr)dr + j;;(3{Xr)a{e)dWr + j!,Jj,J{X,^)Midr,de) 
Yt = g{XT) + f{Qr)dr - Z^dWs - f^^ Ur{e)M{dr, de) 

In fact, there are many methods to develop Malliavin calculus for Levy processes. In our 
paper, we opt for the approach of Sole et al. [25], based on a chaos decomposition in terms 
of multiple stochastic integrals with respect to the random measure M. Adopting notation 
of [12], we will recall the suitable canonical space we adopt to our setting. 

We start by introducing some additional notations and definitions. We assume that the 
probability space (O, P) is the product of two canonical spaces {Q,w x ^/x, Fw x J"^, Wy/ x 
P^) and the filtration F = (-7^t)te[o,T] the canonical filtration completed for P (for details 
concerning this construction, see Section 2 in [12j). 

We consider the finite measure q defined on [0, T] x M by 

q{B)= dt+ e^u{de)dt, B e B{[0,T] x M.). 
Jb{o) Jb' 

where B{0) = {t G [0, T]; (t, 0) G B}, B' = B-B{0) and the random measure Q G [0, T] xR: 

Q{B)= I dWt+ I efl{dt,de), B e B{[0,T] xR). 
Jb{o) Jb' 
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For n € N, a simple function hn = 1eix...x_e„ with pairwise disjoints sets Ei,...,En G 
B{[0,T] X M), we define: 



In{hn)= / h{{ti,ei),...,{tn,en))Q{dti,dei) ■ ... ■ Q{dtn,den). 

j([o,r]xR)" 

We Define the following spaces 

I. Lyg^(M) the space of product measurable deterministic functions h : ([0, T] x M)" 
M satisfying ||/i||?2 < co, where 

T,q,n 

WHh =: / \h{{ti,ei), ...,{tn,en))\'^q{dti,dei) ■ ... ■ q{dtn,den). 



^T,q,n 



([0,T]xR)" 



2. D^'^(]R) denote the space of F-measurable random variables H € L^(M) with the 
representation H = Yl^=Q^n{hn) and satisfying 

oo 

Enn!||/i„||?2 < oo. 

n=0 

3. L^'2(R) denote the space of product measurable and F-adapted processes G : x R — > 
M satisfying 

e(/ |G(s,2/)|2g(ds,dy) I <oo 

yj[o,T]xR y 

G(s, y) G ]D)i'2(E), for 5 - a.e (s, y) G [0, T] x M 
E ( / I y) dz) ) < cx). 

\j([o,T]xK)2 y 

This space is endowed with the norm 

= e(/ |G(s,y)|Vds,dy) ) 

yj[0,T]xIR J 

\j([o,T]xR)2 y 

We should mention that the derivative Dtfi coincide with Dt the classical Malliavin deriva- 
tive with respect to Brownian motion. 

To study the regularity of Z and J7, we shall also introduce the following assumption: 

Assumption 5.1 For 2 < p < | 

1. The generator f has continuous and uniformly bounded first and second order partial 
derivative with respect to x, y and 7. 
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2. For each {x,y,j) e M , dxf{Q),dyf{Q) and 9^/(0) belong to L ' and satisfy 

sup e( [ \Dedif{er)\^dry < oo, (5.2) 
0<6»<T \Je / 

sup sup e( [ \DuD0dif{Gr)dr\A < oo, (5.3) 

where i := x, y, 7. 

There exist a constant K > such that for any e € (M — {0}) , t € [0, T] and 
< 6,u <t <T: 

E\Deg{XT) - DugiXrW < K\e-u\^. (5.4) 

E\De,eg{XT) - DuM^tW < K\e-u\^. (5.5) 

e( [ \Def{er)-Duf{er)A' < K\e-u\i. (5.6) 



\De,ef{er)-Du,ef{Qr)\^y < K\9 - u\2 . (5.7) 
Additionally to Assumption 15.11 we assume that 

Assumption 5.2 For any A > and q > 1, we consider three progressive measurable 
processes {at}o<t<T, {A}o<t<T and {'yt}o<t<T such that: 



Eexp(^xJ^ (\Pr\+lr)dr^ < 



oo, 



sup E \at\'^ + \jt\'^] < oo. 

0<t<T 



Proposition 5.1 Under Assumption \5.^ the discontinuous semi-martingale £t: 

d8t = EtMt + £tlt [ p{e)M{de, dt), {51 

JEe 

has the following properties 



1. E sup < oo, for any n S M. 

0<t<T 

2. The process Zt := satisfies the following linear SDE: 

dZt 



A + it / p\e)m{de) ] dt - -ft / p{e)M {de, dt). 

JE^ J Je^ 

Moreover, we have for any p >2: 

E\Zt-Zs\P <C\t-s\P. (5.9) 
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Proof, f could be written as : 



= exp 







(3r 



[e)m{de) 



dr+ I I 7,.p(e)M((ie, dr) \ . 

JE^ 



Under Assumption 15.21 we get the first assertion. The second assertion is deduced from 
first one, Holder inequality and Burkholder-Davis-Gundy inequality. □ 



The following Theorem constitutes the main tool to prove Theorem [57 

Theoi 

BSDE 



Theorem 5.1 Suppose that 8tXt and arDgXrdr are in M^''^. The following linear 



Yt = g{XT)XT+ I [arXr+ lirYr+-ir'^r]dr - [ ZrdWr 
Jt Jt 



T 



t JE, 



Ur{e)M{dr,de), 0<t<T 



has a unique solution (Y, Z, U) and there is a constant C > such that 

IE|lt -y.r < C\t-s\2 for aUs,te [0,T]. 
Proof. Applying Ito's formula to StYt, we obtain 



(5.10) 



(5.11) 



diStYt) 



-StatXrdt + StZtdWt + £t / {Yatp{e) + Ut{e))M{de, dt). 

'-Be 



Then 



Ft = E {£t,T9{XT)XT + £t,rOlrXrdr/F^ 



where Et^r = Zt£r- 

For < s < t < r, we have: 



-Y,\P < S^'-iE |E {£t,Tg{XT)XT/Ft) - E {£,^Tg{XT)XT / Ts)T 



-3P"^E 



E 



£t,rCirXrdr /Tt^ — ¥, Ss^rOijXrdr / 



(5.12) 



By adapting the argument of Theorem 2.3 in |19j and recall Remark 1 5.1} we can immediately 

p p 

show that Ii < C\t — s| 2 and I2 < C\t — s| 2. □ 



5.1 Malliavin calculus on the Forward SDE 

In this section, we recall well-known properties on forwards SDEs, concerning the Malliavin 
derivatives of the solution of a forward SDE with jump, stated in Nualart [22] in the case 
of SDE without jumps and in Petrou [24] in case of a Levy process. The following theorem 
can be found in 1241. 



20 



Theorem 5.2 

Let X be the solution of forward SDE h5.1\) . Then, for all t € [0,T] and {0,e) G [0,T] x 
(M\{0}), the Malliavin derivatives of X satisfy 



DeXt 



' d^b{Xr)DeXrdr + f [ d^p{Xr)DeXrM{dr,de) 
Je J 

+a{e)l3{Xe)+a{e) [ d^l3{Xr)DeXrdWr, 0<e<t<T. (5.13) 
Je 



And 



Do^eXt = l' De,eb{Xr)dr + [' [ De,eP{Xr)M{dr,de) 

Je J 



+l3{Xe) + a{e) [ De,el^{Xr)dWr, 0<9<t<T. (5.14) 
Je 



For all 9 > t, we have Dq eXt = DgXt = a.s. 



Remark 5.1 Using standard arguments as in the proof of Lemma \4-l\ we can prove the 
following a priori estimate: 



sup E 

o<e<T 



sup DeXt 

0<t<T 



sup E 

0<6»<T,eGM* 



sup De^eXt 

0<t<T 



sup sup E 

O<n<TO<0<T 



sup DuDgXt 

0<t<T 



sup sup E 

0<M<T0<6»<T,eeM* 



sup DuDe,eXt 

0<t<T 



< oo, 



< oo, 



< oo, 



< oo. 



5.2 Malliavin calculus on the Backward SDE 

In this section, we recall some result of Malliavin derivatives applied to BSDE especially 
established in |15] and [12] in the aim to generalize the result of Theorem 2.6 in |19] . 



Theorem 5.3 Assume that assumption \5.1\ hold. There exist a unique solution 
{{Yt, Zt,Ut{e))}Q<t<T,ee{R~{0}) of BSDE ^5J^, such that: 

1. The first version of Malliavin derivative {{DgYf, DgZt, DgUt{e))Q<Q^t<:T^ee{R-{o}) of 
the solution {{Yt, Zt,Ut{e))}Q<t<x,e€{R-{o}) satisfies the following linear BSDE : 

DeYt = d.,g{XT)DeXT+ r f\Xr,YrJr)dr- r DeZrdWs (5.15) 

Jt Jt 

DeUr{e)M{dr, de), 0<e <t<T 



t J El 
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where f{e) := d,f{@)DgXr + dyfiQ)DgYr + ^^/(e)!)^?,. 
Moreover {DtYt)o<t<T is a version of {Zt)o<t<T • 

Zt = DtYt. a.s. (5.16) 

2. The second version of Malliavin derivative {Dg^zYt, Dg^zZt, Dg ,zUt{e))o<e ,t<T,{e,z)&{R-{o}Y 
of the solution (5^t, .^t, ^/t(^;))o<t<T,^G(K-{o}) satisfies the following linear BSDE 

De,,Yt = g{XT + De,zXT)-g{XT) + j\f{er + De,zQr)- f{Qr)]dr (5.17) 
-/ De^zZrdWs- [ I De^zUrie)M{dr,de), 0<e<t<T. 

Moreover (A,e>t)o<i<T,ee(M-{o}) is a version o/ (J/t(2:))o<t<T,ze(lR-{o}) • 

Ut{z) = Dt,eYt a.s. (5.18) 

And for {e,e) G [0,T] x R 

De,zYt = De^zZt = De,zUt{z) = 0, < t < (e,z) G R x (R - {0}). 

3. There exist a constant C > such that for all s,t ^ [0, T]: 



E|Zt-Z^|P < C\t-s\^ (5.19) 
^Tt-TsY < C\t-s\^. (5.20) 

Proof. Existence and uniqueness of solution is similar to Proposition 5.3 in [T2] and 
Theorem 4.1 in [12]. Then we focus our attention to prove inequalities ()5.19p and ()5.20p . 

■ We first prove that E|Zj - < C\t - s\^ . 

Let C > be a constant independent of s and t, whose value vary from line to line. From 
()5.16p we have: 

Zt-Zs = DtYt - DsYs. 

Then 

E\Zt -Zs\P< E\DtYt - D,Yt\P + E\D,Yt - DsY,\p. 
Step 1: Estimate E|A>t - DsYt\P. 
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From Lemma ()6.ip . equation (j5.15p and assumption (j5.4p - ()5.6p . we obtain 



nDtYt-DsYt\P + E / \DtZr-DsZr\'dr 



T \ f 

2, 




+ E{ I j \DtUr{e) - DsUr{e)\^m{de)dr 



< CE(\Dtg{XT) - DsoiXrW 



T 



+CE / \Dtf{r,Xr,Yr,Tr) - Dsf{r,Xr,Yr,rr)\'dr 



< C\t-s\2. 



(5.21) 



step 2: Estimate E\DsYt - DsYs\p 
We recall the expression of £t 

8t = expiyj /3r - ^ 

Denote /J^ = dyf{@r) and 'jr = d^f{Qr)- For any < < t < T, we have: 



dr+ I I 'yrP{e)M{de,dr) 
Je, 



(5.22) 



Wt= £t < pie) d^,f{er)DeXr + d.,yf{er)DeYr + d.,.J{er)DeTr M{de,dr) 

Ue J Ee 

Je J E^ 

f I [dyyfi^r) - p\e)^rdy^f{Qr)] DeYrm{de)dr 

Je J Ee 

f [ [dyyfier) - p2(e)7,a^^/(e,)] I),r,m(de)d^ 
Je J E^ 

From other side, by induction on chain rule: 

De,ef{Q) = f\e + De,ee) - /'(O). (5.23) 
Using this in the previous equality 



+ 



+ 



DgeSt = St exp 



dyf (Qr + De^e^r) " dyf{@r] 



1 



p\e) {d^f{@r + De,eQr)f " 7r rn{de) 



dr 



+ lep{e) + / / p(e) 9^/(6, + De^e^r) - Ir M{de, dr)\-l 
Je J Ee J 

From Proposition 15.11 Assumption 15.11 Holder inequality and Burkholder-Davis-Gundy 
inequality, we can show for any p < q that: 



sup E sup \Dg^eSt\^ < OO. 



(5.24) 
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Now, by Clark-Ocone formula (See Es-Sebaiy and Tudor [16j) on StDqXt- 

+ / / E(Dr,ei£TDeXT)/J'r)M{dr,de) 

JO J 

= E{£TDeXT)+ [ u^^dWr+ [ [ vl,,M{dr,de). 
Jo Jo JSr 



Where 



u^. := Ei^DrSTDeXT + STDrDeXr/Tr 
Thus it remains to prove that 



sup sup \Uj.\^ < oo 



sup sup l^'reT < °0. 



E{Dr,e£TDeXT + £TDr,eDeXT + Dr,e£TDr,eDgXT/Tr ' ^ 



P 

q \ 1 



ee_[0,T] r6[0,T] 

) sup 

ee[o,T] rG[o,r],eei 

By Holder inequality 

< 3"-' f(E|A-,e£Tl")~ {ViDsXtV)' + (E|£rl")~ (l.\D,.„DeXT 

q — p 

+ (E\Dr,e£T\^^ " (E\Dr,eDeXT 

Combining (j5.24p and Remark |5. 11 we deduce that sup^g^-pj sup^g^ ggjg \v^^J^ < oo. Fol- 
lowing the same arguments we conclude that supgg[o,T] suPrG[o,T] Wr\^ < consequence 
£tDqXt belongs to M"^'^. Therefore, by Theorem 15. ll we conclude 

E\D,Yt- DsYs\^ <C\t- s\i. (5.25) 

Now, combining (j5.21|) and (j5.25p . we finally obtain for some constant C > 

E\Zt -Z,f = E\B,Yt - DsYs\^ <C\t- s\-^. (5.26) 

■ We prove that E\Tt - F^l^ <C\t-s\^. 
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By Holder inequality 

Eirf-rjp = E 



/ pie)iUt{e)-Us{e)Mde] 



< \Ut{e)-U,{e)ru{de)j l^J^ \p{e)\—^ u{de) 

< C [ E\Ut{e)-Us{e)\Pu{de) 

= C I E\Dt,eYt-D,,eYsfi^{de) 
Je^ 



< C I E 

IE, 



\Dt,eYt - Ds^eYtf + \Ds,eYt " D^^eYsf y{de). 



step 3.: We prove that E |A,e>^t - Ds,eYtT < C\t - s| 2 . 
Under assumption (jS.Sp . (j5.7p and from Lemma l6.ll 



E\Dt,eYt-Ds,eYt\'' + ^{ 1^ \Dt,eZr - Ds,eZr 



+e( [ [ \Dt,eUr{e) - Ds,eUrie)\^i^{de)d? ' 

\Jt J E^ 

< CE[\Dt,ei-Ds,em 

+CE (^j^ \Dt^J{r,Yr,Ur) - D,J{r,Yr,Ur)\^dr^ ' 

< C\t-s\2. 



(5.27) 



step 4.: We prove that E \Ds,eYt - Ds,eYs\^ < C|i - s| 2 . 
We can write BSDE (|5T7ll as: 



De^.Yt = G{XT)De,eXT + J ae,r[De,eXr + Dg^.Yr + De,,T, 
- [ I De,eUr{e)Midr,de), 

Jt J Ee 



1: 



dr- I Dff^ZrdW, 



(5.28) 



where 



g{XT + De^eXr) - 9{Xt) ^ 

f{@r + De,ee)-f{er) 
De eXr + De eYr + De eTr {^o,.Xr+Do,,Yr+Do,,Vr^O}- 



Then, from Lipschitz continuity of /, we have supo<t<7^ E|a6)^t|^ < 00. It remain to show 



25 



that EtDo^eXt belongs to M^'^. In fact, by Clark-Ocone formula applied to EtDq^i-Xt'- 



+ 



£ E(^Dr,e{£TDe,eXT)/J'r)M{dr, de) 



E(^T-De,e^T) + / uldWr+ [ [ v^^^M {dr , de) , 

Jo Jo J 



with 



Ei^Dr£TDe,eXT + £TDrDe,eXT/T, 
vl^^ := E(^Dr,e£TDe,eXT + £TDr,eDe,eXT + Dr,e£TDr,eDe,eXT/T, 

Following the same argument as Step 2 and using Remark (jS.ip we prove that 



sup E |n^|P + |i}^JP <oo. 
ee[o,T],eeM ^ ' ^ 

Therefore, £tD0 (,Xt belongs to M'^'P. Finally, we apply once again the result of Theorem 
Oto BSDE SUTIElf we get: 



E\Ds,eYt - D,,eYs\ < C\t-S\2. 

The result then follows. 

We now complete the proof of Section 3 

Proof, of Theorem 14.21 We adapt the proof of Theorem 5.2 in [19] 
Let z = n — 0. 

Step 1. : We show that E [supo<i<„ \SZ[f] < C|7r|P-^ 
Denote 

SZl". = Zt. - zi. 



(5.29) 

□ 



Combining (jM]) and (Ii30|) 

E \£t^^Td,,g{X^)DtM - £l^^^^J,g{X^j,)DuX^T / 
E \£t^J\d.,g{X'T) - d^giX^)] Dt,X^ + [a,X|^ - A.^f 



\5Zl\ < 



+ d,giX^)Dt^X^ 



£ti,T - £t.T 
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From Lemma l4.lt inequality (j4.6p and Assumption (A2) 



E sup \6Zl\P < 

0<i<n 



E( sup Dt^X^ 

\0<i<n 



p 
p-l 



PI 



+ 



E 



p -| 1 



E( sup \Dt,X^ - Dt^X^\P 

,0<i<n 



+E 



sup E 

0<i<n 



< Cp 7rP/2+E sup 

0<i<n 



Using the fact that |e^ — e^| < (e^ + e2')|x — y\, leads to 



T 



dr+Tf Mr)p{e)M{de,dr) 



- E / / /3(0p(e)M(rfe,(ir) 



fc=i+l 
n-1 



fc + l r 



1 



/2(r) + - / [/3(r)]V(e)m(de) 



dr 



^+l r 1 

/2(r)-- 
ti L ^ JE^ 



U+i r 

dr+ I / h{r)p{e)M{de,dr) 

ti J Ei: 



It follows that 

h < CE[Dt^X^(^£t^,T + £l^„u) 

[Ur)fp\e)m{de) 

Prom Assumption (A. 3), we have 

\Dt^X^£t,^T\ < \Dt^X^\'exp\ r f2{u)du-l r [ h{u)p'[e)m{de) du 

Uti ^ hi J Et 

+ r f f3{u)pie)Mide,du)] 

Jti J Et ) 

< cf sup \DgX^\] I sup exp\ [ [ h{t)p{e)M{de,du) 

\0<e<T j \0<t<T [Jt JEe 

Similarly, we have: 



< C\ sup \DeX^\ sup exp I / h{t)p{e)M{de,du) 
\G<e<T / Vo<t<T iJt Je^ 
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From other side, for any r > 0, we have by Holder inequahty and Proposition l5.lt 
El sup exp i / / f^{r)p{e)M{de^dr) 

\0<t<T \Jt Je^ 



<E[expl2r [ [ h{r)p{e)M{de,dr) 

I Jo JEe 



xE sup expl-2r I I f-i{r)p{e)M{de,dr) 
\o<t<T I Jo JEe 

< E ( exp I / [r) {e)m{de)dr 



I Jo Jes 



X E sup exp < —2r / / f-i{r)p{e)M{de,dr) 
\o<t<T I Jo Jes 



< oo. 



Thus, for p' E (p. 



E sup If < CE sup \Dt^X^\\ sup exp { / h{t)p{e)M{de,du) 

0<i<n~l \ \0<e<T I \0<t<T [Jt J E^ 



sup / \f2{r)\dr + - sup / 

0<i<n-l Jti ^ 0<j<n-l Hi 



+ sup 

0<i<n-l 



[/sC?^)] P {e)m{de)dr 



ti+l 



h{r)p{e)M{de,dr) 



< C 



E sup |A,^?| 
\ o<e<T 



sup exp\ I I h{t)p{e)M{de,du) 
o<t<T [Jt Je^ 

/ fti+i \ P' 

E sup / \f2{r)\dr 

0<i<n-l \Jti 

+ E sup fr^V [h{r)fp\e)m{de)dr^ 

0<i<n-l \Jt, J Er 



2(p'-p) 



+ E sup 

0<i<n-l 



ti+l 



h{r)p{e)M{de,dr) 



< C h + h + h" 
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We first estimate I3. 

By Holder inequality for r > 1, Jensen inequality and Burkholder-Davis-Gundy inequality: 



< 



< 



E sup 

0<j<ra-l 



I 0<j<n-l 



ti + l 

ti J El 

ti+l 



h{r)p{e)M{de,dr) 



Mr)p{e)M{de,dr) 



rp 



, 0<j<n-l 



[h{r)fp\e)m{de)dr 



V 

rp_ \ ^17 
2 



< CM 



ZtCJ(/T r 

For vr small enough, we take r = — ^7^, then 

if < Clvrf 



p p 

I 2 2!oo J- 



And 



Consequently, 



Step 2. We show that 



E sup \6Zl\P < C\7r\~ """M. 

0<i<n 



E sup \6ri\P < C\7r 

0<i<n 



IP-I 



where ST^ = Tt^ - T^. In fact, 



E sup \6Tl\P = E sup 

0<i<n ' 0<i<n 



p{e)6UlMde] 



< CE sup / pP{e)\6Ul^^\^u{de) 

< CE sup ISUlJ" 



0<i<n 



However, following exactly the same arguments as Step 1, we can prove that: 



p p_ 



E sup |5^7j%|P < C\tt\~ 

0<i<n 



and that 



E P_ 



E sup \5Ttf < C\Tr\~ ^""^ . 

0<i<n 



(5.30) 



(5.31) 
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Step 3.: We show that 



We have 



E sup < C\7r 

0<i<n 



IP-I 



n-1 



\ k=i+l , 

/ n-1 



k=i+l 



Hence by adapted again the argument of [19] to our setting, we have for i = n — 1, n — 2, 0. 

|5y,^| < E \f\®lj - f(&t,jm + \Rl\ + \Sg^X^)\ /Tt] 

\k=i+i J 

where \dg^Xf)\ = \g{X^) - g{X^)\ and 



\Rl 



k=i+l 



For j = n — l,n — 2, i 

\6Y,^\<e(j2 \f{&U-f{&t,+,)m+ sup \Rl\ + \6g-^{X^T)\/j'tA 

Since we know from [19] that KsupQ^^^rp \RJ_\p < C|7r|2, combining this with standard 
estimate of 5X'^ and Lipschitz property of the generator / we have: 

E sup IdY^^'lP 



0<j<n 



< CE 



< CE 



n-1 



E l/(0r,+,)-/(0t.+J|Atfc + sup \Rl\P + \6g^XT) 



\k=i+l 



0<t<T 

\k=i+l J \A;=j+l / \k=i+l J 



+ sup \Rlf^W{XT)f 

0<t<T 



<C{{T-ti)PE sup ISYtZl" + [ \tt\ '^+|7r|^ 

i+l<k<T V 

Using similar recursive methods of Theorem 4.2 in [19] we get estimate: 



E sup Ir^y^'^r < CplTrf 



0<j<n 



(5.32) 
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Finally, combining ()5.30p - (|5.3ip and ()5.32p we close this proof. 



□ 



Proof, of Theorem 14.31 

For any i = n — l,n — 2, 0. Observing that 

max sup E[|y/-y-|2] < C max sup K\\Yf - Yt\' + \Yt - Ytf + \Yt^ - Y,^f 



Combining ()2.10p . ()4.14p with Corollary 2.7 applied to our setting, we have: 
max sup E \\Yf - y.^Pl < C ( a(ef + M ) . 



(5.33) 



Combining (|2.10p and (I5.30p we obtain 
/ VrdRr - V Z^AWt 



E 



< C E 



+ E 





E 

i=0 



VrdRr - I Zf.dWr 




ti+1 



dWr 



< C{a{ef + \TT\ ""^i 



Arguing as above, we obtain 



1 n-^t-^fdt < I E[\Tt-Tl\^ + \Tl-Tlf + \Tl-Tl\^] 



(5.34) 



dt. 



i=0 " 1=0 

From (I2J0]) . ([CT]) and (|OT]) we have: 



V/" E\Tt-Tl\^dt < / E [|r, -r^|2] dr + C7|7r| 

< c([ [ E\Vr-U^.ie)\'^m{de)dr+ [ [ ElVrl"^ m{de)dr 



< C7 c7(e)2 + |7r| 



iog 



(5.35) 



Combining (j5.33p . (|5.34p and (I5.35P we get the required results. 



□ 
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6 Appendix: A priori estimates 



Proof, of Lemma 14.11 

By Jensen inequality 

sup WDeXJ'fP < sup || / d.^b{X^rr)DeX^ndrfP 

0<S<t e<S<t Je r r 

+ S^P-i sup \\a{e) [ d,PiXJ^)DeX^,.dWrfP 
e<s<t Je 

+ 32f-i sup II r [ d,/3{X;„)DeX;„Mide,dr)fP 
e<s<t Je Je^ 

+ 32P-V(e)/3W)I'^ 

Taking expectation in both hand side and using Holder inequality, we obtain: 



E 



sup WDeX^f'P 
e<s<t 



+ E 

< c(e 

< c(e 



j\\DeXl^fPdr 



+ E 



T 



\DeXlA\''Pdr 



\DeXl^fdr 



\DeXU?Pdr 



+ B 



[ sup \\DeXJ„fPdr 

Je 0<u<r 



where B := K [supQ<^<2- ||/3(XJ) Ip^] . Since the constant C doesn't depends on 6 and n, 
we conclude by Gronwall lemma that E [supg<^<7i HDgX^p^] is bounded and therefore 
supo<0<T sup„>i E [sup0<j<r \\DoX]^fP] is finite . 

By the same arguments, we prove that 

sup SUpE sup \\Dg^eX-^ 

o<e<Tn>i e<t<T 



TT |l2p 
t II 



< oo. 



□ 



Lemma 6.1 Let^, G L'?(fi), / : nx[0,T]xRxL'^{E,S,iy; 



measurable, satisfies E |/(t, 0, 0)p(ft < oo and uniformly Lipschitz w.r.t (y, z), such that, 
for some constant K > we have: 

\f{t,yi,ui) - f{t,y2,U2)\ < K{\yi - y2\ + ||tii - n2||), 
for all yi, y2, € M and ui,U2 € L^{E, £, v; M). 

Then, there exist a unique triple (y, Z, U) e B'^ solution to BSDEs i f73|) . Moreover, For 
q >2, we have the following a priori estimate: 



E 



sup \Yt\i 

0<t<T 



+ E| y \Zt\ 



+ E( 1^ [ \Utie)\'^u{de)dt 
lo Je 



T \ 2 ' 

< C(E|er + E( |/(t,0,0)|2dt 



(6.1) 
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Proof. Existence and uniqueness of solution of BSDEs with jump are proved in [3] and 
the estimate (j6.ip is a direct consequence of proposition 2.2 in the same reference, with 
(/',Q') = (0,0). □ 
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